Prospects for Constraining Neutrino Mass Using Planck and Lyman- Alpha Forest 
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In this paper we investigate how well Planck and Lyman- Alpha forest data will be able to constrain 
the sum of the neutrino masses, and thus, in conjunction with flavour oscillation experiments, be 
able to determine the absolute masses of the neutrinos. It seems possible that Planck, together with 
a Lyman-Q survey, will be able to put pressure on an inverted hierarchial model for the neutrino 
masses. However, even for optimistic assumptions of the precision of future Lyman-Q datasets, it 
will not be possible to confirm a minimal-mass normal hierarchy. 



I. INTRODUCTION 



The determination of absolute neutrino masses is a key 
scientific goal for the coming decade. Neutrino flavour os- 
cillation detections have shown that neutrinos do indeed 
have mass but unfortunately cannot determine their ab- 
solute masses. Particle physics experiments (e.g. tritium 
beta decay or neutrinoless double beta decay) offer the 
most direct probe of neutrino masses, but reaching limits 
of less than 1 eV is formidably challenging Sub-eV 
neutrino masses can also be probed indirectly via their 
effects on the energy density of the Universe and large 
scale structure. For this reason there has been consid- 
erable interest on constraints on neutrino masses from 
various current and future cosmological probes (see [2| 
for an excellent and comprehensive review). 

In this paper we perform a detailed Markov Chain 
Monte Carlo analysis to assess the sensitivity of the 
Planck satellite [Sj] to neutrino masses. This problem 
is topical because Planck is scheduled for launch in late 
2008 and should provide all-sky maps of the cosmic mi- 
crowave background (CMB) of unprecedented precision. 
Planck will operate at nine frequencies from 30 to 857 
GHz, seven of which will be sensitive to polarization. The 
sensitivity and frequency coverage of Planck should allow 
accurate subtraction of foregrounds resulting in an essen- 
tially complete reconstruction of the CMB temperature 
signal over a large area of the sky together with signifi- 
cant new information on its polarization. For a summary 
of the Planck instruments and its scientific programme, 
including constraints on neutrino masses, see 0]. For as- 
sessments of Planck constraints on neutrino masses from 
different perspectives and using a variety of complemen- 
tary astrophysical data see [1, [1, H, 0, Hi ■ 

Neutrino masses affect both the cosmic history and 
structure formation. One main effect is a reduction of 
power below a characteristic wavenumber corresponding 
to the Hubble scale when the neutrinos become non- 
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(e.g. [l,[i3j[iH)- Planck will not directly probe the damp- 
ing of fluctuations caused by massive eV-scale neutrinos 
(except possibly via weak lensing of the CMB). Never- 
theless Planck will be vital since it will break parameter 
degeneracies that would otherwise exist in other cosmo- 
logical datasets. 

Observations of Lyman-a absorption in quasar spec- 
tra constrain the amplitude and slope of the cosmological 
matter power spectrum at redshifts between two and four 
on comoving megaparsec scales. Compared to galaxy red- 
shift surveys Lyman-a surveys probe the matter power 
spectrum closer to the linear regime and sample smaller 
scales, giving a long lever arm when combined with ob- 
servations of the CMB. Furthermore, the galaxy power 
spectrum is difficult to relate to the matter power spec- 
trum on small scales where the fluctuations are highly 
non-linear and the damping effects of eV-mass neutrinos 
are most significant. Hence it is natural to consider com- 
bining Planck with a Lyman-a survey in an effort to as- 
sess what cosmology might contribute to to the determi- 
nation of neutrino masses over the next decade. Already, 
the tightest cosmological limits on neutrino masses come 
from combining Lyrnan-a datasets with CMB and other 
datasets (e.g. [13, IS)- principle an ultra-large galaxy 
redshift survey, as anticipated from the Square Kilometer 
Array (SKA), might be able to measure the small effects 
of sub-eV scale neutrinos on the galaxy power spectrum 
on scales > 0.02/iMpc~^ (see Q). This is discussed fur- 
ther in Sections IIXI and 1x1 

We first consider constraints on neutrino masses from 
the current Sloan Digital Sky Survey Lyman-a data (de- 
noted SDSSLYA below). We then consider the improve- 
ment in these contraints arising from a hypothetical 
percent-level determination of the power spectrum from 
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a future Lyman-a survey. Neutrino flavour oscillation ex- 
periments are consistent with one of two minimal values 
for the sum of neutrino masses, either 0.056 eV or 0.095 
eV. The former value comes from assuming the neutri- 
nos fall into a normal hierarchy, with the mass of the 
intermediate-mass neutrino closer to that of the lightest 
neutrino than that of the heaviest one. The latter value 
comes from assuming the neutrinos fall into an inverted 
hierarchy, with the mass of the intermediate-mass neu- 
trino closer to that of the heaviest neutrino than that of 
the hghtest one (see e.g. For a total mass signifi- 

cantly greater than these values, the mass splittings are 
much smaller than the individual masses and the neutri- 
nos are said to be degenerate. These observations sug- 
gests two clear goals for a cosmological determination of 
neutrino mass: firstly to see whether neutrino masses are 
degenerate or not and secondly to differentiate between a 
normal and an inverted hierarchy. In addition, it may be 
possible to put useful constraints on physics beyond the 
standard model that impinges on neutrino mass; e.g. the 
study of thermal leptogenesis in that gives an upper 
bound on the mass of a light neutrino of 0.1 eV. 



II. COSMOLOGICAL IMPLICATIONS OF 
NEUTRINO MASS 

The evolution of both the background universe and 
perturbations within it are sensitive to the neutrino mass 
spectrum. For reviews see 0, [l^. One main effect of 
massive neutrinos is a suppression of power on small 
scales (Eq. ([1])) that is roughly proportional to the neu- 
trino fraction of the matter content of the universe. How- 
ever the effect can be partially degenerate with changes in 
other cosmological parameters, so for robust results un- 
certainties in other parameters should be marginalized 
out. 



III. MASSIVE NEUTRINOS IN CAMB AND 
COSMOMC 

Our analysis makes use of the software packages 
CAMB^ [ill and CosmoMC^ CAMB calculates the 

linear-theory CMB power spectrum and optionally multi- 
ple matter power spectra for a given cosmological model. 
CosmoMC uses Markov Chain Monte Carlo to sample 
from the posterior distribution of cosmological param- 
eters from a given likelihood function. The likelihood 
function uses theoretical calculations from CAMB in 
combination with real (or mock) datasets. 

The standard set of six parameters commonly used 
is {rib/i^, f^dm/i^, T, rig, ln(10^° As)}, where fib is the 
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baryon density divided by the critical density, fldm is the 
dark matter density (including potential massive neutri- 
nos) divided by the critical density, h is the Hubble con- 
stant in units of 100 km s~^ Mpc~^, 9 is the acoustic 
horizon angular scale, t is the optical depth from reion- 
ization, Ag and rtg are the amplitude and spectral index 
(at a fiducial wavenumber of 0.05 Mpc~^) of the pri- 
mordial adiabatic scalar curvature perturbation power 
spectrum V{k). We follow this usage here, and use the 
shorthand {6} to denote these six parameters. In this 
paper, the dark energy is assumed to be a cosmological 
constant. Massive neutrinos are introduced via the pa- 
rameter fi, = rii^/ridm, where J7y is the massive neutrino 
density divided by the critical density today. The sum 
^ rrii, of the neutrino masses is related to fi, by [2] 

^TO^ « 93.12r2dm/iV^ eV. (2) 

Throughout this paper we assume only the usual three 
neutrino species. We either take them all to have the 
same mass ^mj//3, or take two of them to be mass- 
less and one to be massive with mass w (= ^ in 
this case). In the first case we denote the standard six 
parameters along with ^ rrii, by {6 -I- S} and in the sec- 
ond case we write {6 -I- m}. We shall also sometimes 
consider a running scalar spectral index, with running 
rij-un = {d/dhik)^hiV{k) and assumed constant in k. 

The publically available version of CAMB has recently 
been upgraded by one of us (AL) to handle arbitrary mass 
splittings. Fig. [T] shows how the matter power spectrum 
is sensitive to assumptions about the neutrino masses. 
For this work we have also modified CosmoMC to allow 
two neutrinos to be massless and one to be massive. In 
all cases CosmoMC uses fi, as its base parameter, with m 
or ^ as derived parameters. 

IV. WMAP3 AND MASSIVE NEUTRINOS 

In Ref. [13] the WMAP team present constraints on neu- 
trino properties based on the 3-year WMAP data [H, [l^ 
either with or without other datasets. Using WMAP data 
alone, they find '^mu < 1.8 eV at 95% confidence. 
Along with either SDSS or 2dFGRS galaxy redshift data, 
they find X^'^i' < 1-3 eV or 0.9 eV at 95% confidence 
respectively. We ran seven-parameter {6 -I- S} chains 
against wmap and found 'Y^rrii, < 1.7 eV at 95% confi- 
dence (using version 2 of the wmap team likelihood code) 
in good agreement with [12| . 

V. WMAP3 AND SDSSLYA 

A recent paper obtains an impressive 95% upper bound 
for of oiily 0-1''' sV, using a combination of CMB, 

galaxy and Lyman-a data [l3l | . Their routine to calculate 
the likelihood for a model in light of the SDSSLYA data 
has been made publicly available at [2^. The CMB and 
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FIG. 1: A plot illustrating the dependence of the mat- 
ter power spectrum on the neutrino mass spectrum. The 
power spectra of three models are illustrated relative to the 
power spectrum for three equal-mass neutrinos with ^2, ~ 
0.059 eV. The upper solid blue curve corresponds to mass- 
less neutrinos, the dotted curve corresponds to the mini- 
mal normal hierarchy, and the lower curve corresponds to 
the minimal inverted hierarchy. Also shown is a model with 
^rrii, = 0.13 eV and other cosmological parameters changed 
so that the model is nearly degenerate in likelihood with the 
fiducial one against future Planck and Lyman-Q datasets. 

other data are somewhat in tension, the small-scale data 
favouring a higher overall normalization of the power 
spectrum than the CMB. Since adding neutrino mass only 
lowers the power on small scales, the combined datasets 
prefer no neutrino mass at all. The tension between the 
datasets could be merely a statistical fluctuation: the 
two datasets happen to give a tighter mass constraint 
than expected from most a priori possible realizations 
of the data. Alternatively it could indicate that there 
is some inconsistency in the datasets, for example one, 
or both, having some unaccounted-for systematic error. 
Or it could indicate that our modelling is too simplified, 
e.g. the primordial power spectrum varies strongly with 
wavenumber. 

We performed a seven-parameter {6 -I- E} joint anal- 
ysis of the WMAP and SDSSLYA data, and obtained a 
95% upper bound on J2 "^i^ of 0.39 eV (v2 wmap code) 
or 0.35 eV (vl wmap code), intermediate between that 
of [l^ and those of the wmap team using galaxy sur- 
vey data mentioned above. Note that we used the same 
WMAP and Lyman-a forest data and code as How- 
ever we did not use the additional CMB galaxy survey 
and supernovae data used by Ref. [13] ; including these 
datasets gives them a significantly tighter constraint than 
from just Lyman-a and WMAP alone. 

To investigate the possible inconsistency of the wmap 
and SDSSLYA data, we constructed a fake WMAP dataset, 
keeping the noise the same as for the real data but replac- 
ing the Ci 's themselves by those of the 6-parameter model 
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FIG. 2: A 2-D marginalized likelihood contour plot indicat- 
ing the possible inconsistency of the wmap and SDSSLYA data. 
68% and 95% confidence intervals are illustrated for the fol- 
lowing four datasets (from broadest to tightest): blue, WMAP 
data alone; yellow, "faked" wmap data alone; green, SDSSLYA 
and wmap data; red, sdsslya and "faked" wmap data. 



that is the best fit to the wmap and SDSSLYA datasets 
taken together. With {6 + E} chains, faked wmap alone 
gives a 95% upper limit to ^ of 1.5 eV, compara- 
ble to that from real wmap. However, faked wmap with 
SDSSLYA gives a 95% upper limit to J^^i' of O-^O eV, 
double that of real wmap with SDSSLYA; see Fig. [51 A 
possible interpretation of this result is that the WMAP 
+ SDSSLYA limit on the neutrino mass is spuriously low, 
by a factor of about two, because of some unidentified 
systematic error in one or both datasets. However, as 
mentioned above, the tension between the two datasets 
seen in Fig.[2]could be a statistical fluke or an inadequacy 
in the theoretical model. 



VI. MOCK PLANCK DATA FOR COSMOMC 

In this section we describe how we construct very sim- 
plified mock Planck data for CosmoMC, used for fore- 
casting future constraints on the neutrino masses. We 
consider only the linear-theory CMB power spectrum. If 
the CMB lensing signal (via the power-spectrum of the 
weak lensing deflection field) can also be used, signifi- 
cantly better constraints might be obtainable than by 
using the CMB power spectrum alone. This is discussed 
by Lesgourgues et al. in Q, who find that including the 
weak lensing deflection field from an idealised Planck ex- 
periment (unlikely to be realised in practice) leads to lim- 
its on the neutrino masses comparable to those discussed 
below from combining Planck with Lya data. 

First of all we assume an underlying cosmology, tak- 
ing parameters from the best-fit six-parameter models 
coming from either wmap alone or wmap and SDSSLYA. 
Next we run the CAMB software to generate a theoreti- 
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cal "prior" power spectrum distribution with mean Cf^ 
based on the input cosmological parameters. For com- 
parison purposes we sometimes include one non-zero neu- 
trino mass 0.06 eV when generating this power spectrum. 

The log-Ukehhood for some model with power spec- 
trum Ci averaged over sky realizations turns out, up 
to an irrelevant constant, just to be the log-likelihood 
of those Cj's evaluated taking the sky power spectrum 
to be its ensemble average Cf^. Hence we do not need 
to make a specific realization of the sky or numerically 
average over many of them for our forecasting; we just 
imagine the data to have exactly the ensemble average 
power spectrum. This procedure gives error bars consis- 
tent with those obtained from most actual realizations, 
but has the advantage that the maximum likelihood pa- 
rameters should be at their true values rather than mov- 
ing around between different realizations. (See [2l|, [g^l 
for related discussions.) 

We then create the mock dataset for use with 
CosmoMC [23]. The sky power spectrum Ci is set to 
be the prior spectrum described in the previous para- 
graph. We work on the full sky and convolve with a 
beam window function and add noise (see e.g. [l^] for 
the procedure). The beam is assumed to be Gaussian 
and symmetric, and the noise is assumed to be white and 
uniform across the sky (we neglect any foregrounds) . The 
beam width and pixel noise are chosen to approximate 
those of the Planck 143 GHz High Frequency Instrument 
channel [11 . 

We label the prior power spectra derived from current 
data as follows: from the six-parameter fit to wmap 
alone assuming massless neutrinos, from the six- 

parameter fit to WMAP and SDSSLYA assuming massless 
neutrinos, Cggg from the six-parameter fit to wmap as- 
suming one neutrino mass of 0.06 eV and Cp^g from the 
six-parameter fit to wmap and SDSSLYA with one neu- 
trino mass of 0.06 eV. 



VII. PLANCK ALONE 

Running CosmoMC for seven-parameter chains {6 + S} 
against the no-neutrino-mass, WMAP-only model 
yields ^ O-'^l at 95% confidence. Thus Planck 

by itself should do at least twice as well as wmap in con- 
straining the sum of neutrino masses. Of course we have 
assumed an idealized beam and noise model for Planck^, 
but on the other hand we have assumed the sensitiv- 
ity for only one frequency channel. The precision from 
CMB data alone (without CMB lensing) is limited since 
neutrinos light enough to be effectively massless at re- 
combination have almost no effect on the CMB power 
spectrum. 



* Residual 'striping' noise for Planck should be ignorable in com- 
parison to white noise, see [251 . 



We also ran eight-parameter {6 + 'E + n^^n} chains with 
both neutrino mass ^ and scalar spectral index run- 
ning TT-run. Despitc thoughts that both might affect the 
power spectrum on small scales in a qualitatively similar 
way, we found that these two parameters were not in fact 
degenerate and that the 95% upper limit on ^ rrii, was 
only moderately weakened to 0.87 eV. 

VIII. PLANCK AND SDSSLYA 

We now consider how well Planck might do in con- 
junction with SDSSLYA data. Running {6 -I- S} chains 
against and SDSSLYA, we obtain the tight constraint 
< 0.10 eV at 95% confidence, three times tighter 
than WMAP and SDSSLYA suggesting that Planck data 
together with existing Lyman-a data may be capable of 
placing severe pressure on an inverted hierarchy. The 
same constraint on ^ is obtained using Cg^g in place 
ofCW. 

However, there is a concern that this constraint might 
be artificially tight simply because of the possible dis- 
crepancy between the WMAP and SDSSLYA datasets dis- 
cussed in Section |Vl To address this, we ran {6 -I- E} 
chains against and SDSSLYA, rather than and 

SDSSLYA. Although the SDSSLYA data effectively enters 
twice, this procedure should give an indication of what 
might happen if the tension between the wmap and 
SDSSLYA data is caused by systematic errors. This yields 
^ m,y < 0.27 eV, a significantly weaker constraint. Thus 
we see that that bound on the neutrino mass is highly 
sensitive to the assumed input model. The ability of 
Planck and SDSSLYA to constrain an inverted hierarchy 
therefore depends on which of these input models is closer 
to the truth. 



IX. PLANCK AND PERCENT LEVEL 
MEASURES OF THE MATTER POWER 
SPECTRUM 

In this Section we analyse what might be learned from 
Planck and new Lyman-a surveys of g reater statistical 
power than SDSSLYA. For example, [2^ investigates how 
an extended Lya survey might perform in constraining 
dark energy and curvature, assuming an experimental 
configuration with characteristics similar to that pro- 
posed for galaxy baryonic acoustic oscillation surveys. 
Here we take a very simple approach and consider a sur- 
vey that would be able to measure the matter power spec- 
trum at one or more effective redshifts and at one or more 
scales to better than five percent accuracy. 

Since Lyman-a surveys effectively measure distances in 
velocity units (see e.g. [23|), we choose our scales likewise. 
We consider the following hypothetical datasets: 

1. Pi®^, consisting of a single data point at an effec- 
tive redshift of 3 and wavenumber 0.009 s/km, with 
a 1% fractional error. 
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2. Pi%^ , consisting of three data points at an effec- 
tive redshift of 3 at wavenumbers 0.002 s/km, 0.009 
s/km and 0.02 s/km, with 1% fractional errors, 

3. Pi%^ , consisting of three data points at each of 
three redshifts of 2, 3 and 3.5 at wavenumbers of 
0.002 s/km, 0.009 s/km and 0.02 s/km, with 1% 
fractional errors, and 

4. P^%^ , consisting of three data points at each of 
three redshifts of 2, 3 and 3.5 at wavenumbers of 
0.002 s/km, 0.009 s/km and 0.02 s/km, with 5% 
fractional errors. 

The datasets are constructed by evolving the z = 
matter power spectrum output from CAMB at the appro- 
priate scale back to the appropriate redshift using the 
standard formula for the growth of linear inhomogeneities 
with the appropriate parameters for the assumed back- 
ground model. This "data" is fed into a version of the 
lya.f90 module of CosmoMC in order to perform the 
likelihood calculation for models. 

Running {G-t-m-t-nrunj-parameter chains against C^qq 
and faked Lyman-a datasets yields the limits shown in 
Table |T] for ^ at 95% confidence. 

Corresponding one-dimensional likelihood plots are 
shown in Fig. [31 Note that all curves peak in the vicin- 
ity of the added neutrino mass of 0.06 eV. Fig. |4] shows 
how adding future Lyman-a data to Planck data breaks 
degeneracies and thus substantially improves the limits 
shown in Fig. [21 However, none of the curves in Fig. [31 
tend to zero as the neutrino mass tends to zero and thus 
none of the dataset combinations is capable of unam- 
biguously detecting neutrino mass in the minimal normal 
hierarchy model. 

In this paper we have focussed on combining CMB and 
Lyman-a data for the reasons given in the introduction. 
In light of the above results we also considered the ad- 
ditional degeneracy-breaking effects that a future galaxy 
survey might provide. Combined with our most opti- 
mistic Pi%^ Lyman-a dataset along with Planck, such 
a galaxy survey would have to effectively measure ^Imh 
to better than 2% in order to yield a 95% confidence 
positive detection of neutrino mass for the minimal hi- 
erarchy (assuming the improved constraint comes from 
degeneracy-breaking alone). For comparison, the SDSS 
and 2dF galaxy surveys constrain Qmh to an accuracy 
of about 10% [28), i29,] and so substantially larger red- 
shift surveys would be required to constrain the shape of 
the matter power spectrum to the level required to con- 
strain a minimal hierarchy. At present, the best prospect 
seems to be a large-scale galaxy survey of ~ 10^ galaxies 
detected with the SKA over the redshift range — 1.5 [sj. 

X. CONCLUSIONS 

In this paper we have studied how Planck, in com- 
bination with a Lyman-a based measure of power on 



TABLE I: A Table showing the limits on rn,^ obtained using 
the assumed futuristic Lyman-a datasets (denoted by the P's) 
and Planck dataset (denoted by the C). 
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FIG. 3: A plot of the marginalized likelihoods for a single 
neutrino of mass m with assumed future datasets as discussed 
in the text. AH curves use the Co^g Planck dataset. As for 
the Lyman-a dataset used, black (solid) corresponds to Pi%^ , 
red (dot-dash) to Pl%^ , green (short-dash) to Pf%^ and blue 
(long-dash) to Pf|^ 
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FIG. 4: A 2D contour plot indicating how a partial parame- 
ter degeneracy using only Planck data is lifted when Lyman- 
a data is added. 68% and 95% confidence intervals are il- 
lustrated for the following three datasets (from broadest to 
tightest): blue, Planck alone; green, Planck with Pl^J'- red, 
Planck with P'fS^ . 
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megaparsec scales, might perform in constraining neu- 
trino masses. 

We find tliat Planck, in combination witli existing 
Lyman-a data, should be able to put significant pressure 
on the inverted hierarchy model. Some of the allowed 
parameter space for thermal leptogenesis models should 
also be constrained. These limits can be tightened by 
using more powerful Lyman-a data. However, even for 
the ambitious Lyman-a datasets assumed in Section HXl 
we conclude that it is unlikely that Planck will be able 
to positively detect a minimal-mass normal hierarchy. 

Other cosmological data, in particular lensing of the 
the CMB, may improve the neutrino mass constraints, 
though extracting an accurate lensing signal from real- 
istic Planck data is likely to be challenging and needs 
further investigation. For idealized Planck data, [7| con- 
clude that the weak lensing deflection field can improve 
the neutrino mass limits from Planck alone to a la limit 
of 0.13 eV, comparable to our forecasts for Planck com- 
bined with SDSSLYA. In principal, a sufficiently large 
galaxy redshift survey, such as envisaged for the SKA, 
in combination with Planck could probe a minimal-mass 
normal hierarchy. This has been considered in some de- 
tail in [8^. Apart from the long timescale involved for 
such a survey (probably well beyond the next decade) it 
may prove difficult to relate the galaxy power spectrum 
to the underlying matter power spectrum to the required 
accuracy. (See [29[ for empirical evidence of scale depen- 
dent bias in the galaxy distribution over the wavenumber 
range 0.01 < k < 0.15/iMpc~^.) Whether the precision 
envisaged by [1] can be achieved remains to be seen. 



The absolute values of the neutrino masses would of- 
fer important insights into physics beyond the standard 
model. There is a widespread hope that cosmological 
probes will be able to constrain neutrino masses to a 
precision better than the normal hierarchy characteris- 
tic mass of 0.06 eV. However, the detailed calculations 
presented here suggest that we should be more sanguine. 
The cosmological detection of 0.06 eV neutrinos would 
require extremely large cosmological datasets, free of sys- 
tematic errors, in addition to Planck. Furthermore, the 
cosmological limits are dependent on physical assump- 
tions (e.g. featureless varying power spectrum, and fixed 
dark energy) that may be difficult to justify experimen- 
tally. A convincing detection of a neutrinos mass < 0.1 
eV will require, at the very least, consistency between a 
number of independent cosmological datasets. 
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